clear all

use "C:\Users\USUARIO\Desktop\Tesis\Tema de I+D\Base Stata\DoFile y Base (avances)\Base y do file 4\base_despues_de_test_validacion.dta", clear

**errores en las variables - redefinidas adelante

drop PHD
drop maestria
drop especialista
drop tercer_nivel
drop tecnico
drop secundaria
drop primaria


*ajuste 1 - definiendo 2 tipos de cooperantes - coopeantes de proceso y conocimiento 

gen coop_proc=0
replace coop_proc=1 if col_clientes==1
replace coop_proc=1 if col_compet==1
replace coop_proc=1 if col_prove==1
replace coop_proc=1 if col_emp_rel==1
replace coop_proc=1 if col_emp_dg==1

label var coop_proc "cooperantes de proceso (clientes-compet-pove-empr_relac_empr_del_grupo)"

tab coop_proc

*ajuste 2 - proporcion cooperantes de proceso / cooperantes de conocimiento 

gen tot_coop_proc=col_clientes+col_compet+col_prove+col_emp_rel+col_emp_dg
gen tot_coop_conoc=col_consul + col_univ + col_lab + col_org_pub + col_ofi_pi

label var tot_coop_proc "total de agentes del tipo de proceso con el que cooperan"
label var tot_coop_conoc "total de agentes del tipo de conocimiento con el que cooperan"

tab tot_coop_proc
tab tot_coop_conoc


replace tot_coop_proc = 0.01 if tot_coop_proc==0
replace tot_coop_conoc = 0.01 if tot_coop_conoc==0

gen prop_proc_conoc=tot_coop_proc/tot_coop_conoc

*problema porque se divide por cero - CORREGIDO PONIENDO 0.00001

label var prop_proc_conoc "proporci�n de cooperantes de proceso / cooperantes de conocimiento"



*Uniendo Peque�a y Microempresa

gen Peq=0
replace Peq=1 if micro==1
replace Peq=1 if peque�a==1

label var Peq "1 si es Peque�a o Microempresa"



*EXPERIENCIA 

gen exp=.
replace exp=2014-ii_6 if a�o == 2014
replace exp=2011-II_6 if a�o == 2011

gen exp_2=exp*exp

label var exp "a�os de vida de la empresa"
label var exp_2 "a�os de vida de la empresa al cuadrado"



*DESTINO DE LOS PRODUCTOS 

gen ecuador=0
replace ecuador=1 if II_10_a==1 
replace ecuador=1 if ii_10_a==1

gen latinoamerica=0
replace latinoamerica=1 if II_10_b==2
replace latinoamerica=1 if ii_10_b==2

gen eeuu_can=0
replace eeuu_can=1 if II_10_c==3
replace eeuu_can=1 if ii_10_c==3

gen europa=0
replace europa=1 if II_10_d==4
replace europa=1 if ii_10_d==4

gen asia=0
replace asia=1 if II_10_e==5
replace asia=1 if ii_10_e==5

gen africa=0
replace africa=1 if II_10_f==6
replace africa=1 if ii_10_f==6

gen oceania=0
replace oceania=1 if II_10_g==7
replace oceania=1 if ii_10_g==7

label var ecuador "1 los productos estuvieron dirigidos a Ecuador "
label var latinoamerica "1 los productos estuvieron dirigidos a Latinoamerica "
label var eeuu_can "1 los productos estuvieron dirigidos a EEUU y Canada  "
label var europa "1 los productos estuvieron dirigidos a Europa  "
label var asia "1 los productos estuvieron dirigidos a Asia "
label var africa "1 los productos estuvieron dirigidos a Africa "
label var oceania "1 los productos estuvieron dirigidos a Oceania  "

table ecuador, contents(freq ) by(a�o)
tab II_10_a
tab ii_10_a
 
table latinoamerica, contents(freq ) by(a�o)
tab II_10_b
tab ii_10_b
 
table eeuu_can, contents(freq ) by(a�o)
tab II_10_c
tab ii_10_c

table europa, contents(freq ) by(a�o)
tab II_10_d
tab ii_10_d

table asia, contents(freq ) by(a�o)
tab II_10_e
tab ii_10_e

table africa, contents(freq ) by(a�o)
tab II_10_f
tab ii_10_f

table oceania, contents(freq ) by(a�o)
tab II_10_g
tab ii_10_g


*Logarirmo de total exportaciones y gastos en I+D)

gen l_ga_IyD = log( ga_IyD+1)
gen l_tot_exp=log( tot_exp+1)


*Clasificacion sectores seguin ciiu a 1 digito

gen sect_manuf=0
replace sect_manuf=1 if ciiu1=="C" 

gen sect_minas=0
replace sect_minas=1 if ciiu1=="B"

gen sect_com_bie_serv=0
replace sect_com_bie_serv=1 if ciiu1=="A"
replace sect_com_bie_serv=1 if ciiu1=="D"
replace sect_com_bie_serv=1 if ciiu1=="E"
replace sect_com_bie_serv=1 if ciiu1=="F"
replace sect_com_bie_serv=1 if ciiu1=="G"
replace sect_com_bie_serv=1 if ciiu1=="H"
replace sect_com_bie_serv=1 if ciiu1=="I"
replace sect_com_bie_serv=1 if ciiu1=="J"
replace sect_com_bie_serv=1 if ciiu1=="K"
replace sect_com_bie_serv=1 if ciiu1=="L"
replace sect_com_bie_serv=1 if ciiu1=="M"
replace sect_com_bie_serv=1 if ciiu1=="N"
replace sect_com_bie_serv=1 if ciiu1=="O"
replace sect_com_bie_serv=1 if ciiu1=="P"
replace sect_com_bie_serv=1 if ciiu1=="Q"
replace sect_com_bie_serv=1 if ciiu1=="R"
replace sect_com_bie_serv=1 if ciiu1=="S"

label var sect_manuf "1 si pertenece al sector manufacturero "
label var sect_minas "1 si pertenece al sector minas "
label var sect_com_bie_serv "1 si pertenece al sector comercio, bienes y sevicios "

			*******************ajustes 16/06/2021
			
*Publica o Privada

gen publica=0

replace publica=1 if ii_3==1
replace publica=1 if II_3==1

label var publica "1 si la empresa es publica"


*Obst�culo - Falta de personal calificado en la empresa

gen obs_pers_cal=0
replace obs_pers_cal=1 if X_1_3_a == 1  
*replace obs_pers_cal=1 if X_1_3_a == 2
replace obs_pers_cal=1 if x_1_3_a == 1
*replace obs_pers_cal=1 if x_1_3_a == 2

label var obs_pers_cal "1 si la importancia a la falta de pers. calif. en la empresa es alta "

*Obst�culo - Dificultad de encontrar socios

gen obst_socio=0

replace obst_socio = 1 if X_1_3_e == 1
*replace obst_socio = 1 if X_1_3_e == 2
replace obst_socio = 1 if x_1_3_e == 1
*replace obst_socio = 1 if x_1_3_e == 2

label var obst_socio "1 si la importancia a la dificultad de encontrar socios es alta "

*Obst�culo - Dificultad COSTOS ALTOS

gen obst_cost_alt=0

replace obst_cost_alt = 1 if x_1_1_c == 1
replace obst_cost_alt = 1 if X_1_1_c == 1

label var obst_cost_alt "1 si la importancia a la dificultad por costos altos es alta "

*Obst�culo - Dificultad Fondos en la empresa

gen obst_fondos=0

replace obst_fondos = 1 if x_1_1_a == 1
replace obst_fondos = 1 if X_1_1_a == 1

label var obst_fondos "1 si la importancia a la dificultad por fondos en la empresa es alta "


*Obst�culo - Dificultad Financiamiento externo alto

gen obst_fin_ext=0

replace obst_fin_ext = 1 if x_1_1_b == 1
replace obst_fin_ext = 1 if X_1_1_b == 1

label var obst_fin_ext "1 si la importancia a la dificultad por costos altos de financiamiento externo es alta "


*Obst�culo - Dificultad Falta de personal en el pa�s

gen obst_pers_pais=0

replace obst_pers_pais = 1 if x_1_3_b == 1
replace obst_pers_pais = 1 if X_1_3_b == 1

label var obst_pers_pais "1 si la importancia a la dificultad de tener personal calificado en el pa�s es alta "


*Obst�culo - Informaci�n sobre tecnolog�a

gen obst_inf_tecno=0

replace obst_inf_tecno = 1 if x_1_3_c == 1
replace obst_inf_tecno = 1 if X_1_3_c == 1

label var obst_inf_tecno "1 si la importancia a la dificultad de no tener informacion de tecnologia es alta "

*Obst�culo - Dificultas encontrar socios

gen obst_enc_socio=0

replace obst_enc_socio = 1 if x_1_3_e == 1
replace obst_enc_socio = 1 if X_1_3_e == 1

label var obst_enc_socio "1 si la importancia a la dificultad de encontrar socios es alta "

*Porcentaje de financiamiento de "Apoyos gubernamentales"
gen porc_apoy_gob=.
replace porc_apoy_gob = vi_1_a if a�o==2014
replace porc_apoy_gob = VI_1_a if a�o==2011


*Porcentaje de financiamiento de "Banca privada"
gen por_apoy_bancapr=.
replace por_apoy_bancapr = vi_1_b if a�o==2014
replace por_apoy_bancapr = VI_1_a if a�o==2011

*Porcentaje de financiamiento de "Recursos provenientes del exterior"
gen por_apoy_exter=.
replace por_apoy_exter = vi_1_c if a�o==2014
replace por_apoy_exter = VI_1_a if a�o==2011

*Porcentaje de financiamiento de "Recursos propios"
gen por_apoy_rec_prop=.
replace por_apoy_rec_prop = vi_1_d if a�o==2014
replace por_apoy_rec_prop = VI_1_a if a�o==2011

*Porcentaje de financiamiento de "Otras"
gen por_apoy_otros=.
replace por_apoy_otros = vi_1_e if a�o==2014
replace por_apoy_otros = VI_1_a if a�o==2011


label var porc_apoy_gob "Porcentaje de financiamiento de Apoyos gubernamentales"
label var por_apoy_bancapr "Porcentaje de financiamiento de Banca Privada "
label var por_apoy_exter "Porcentaje de financiamiento de Recursos del Exterior "
label var por_apoy_rec_prop "Porcentaje de financiamiento de Recursos Propios "
label var por_apoy_otros "Porcentaje de financiamiento de Otros "


************************termina ajuste 16/6/2021



*Objetivo - Mejor� la calidad de bienes o servicios

gen obj_cal_bie=0

replace obj_cal_bie = 1 if XIII_1_e == 1
replace obj_cal_bie = 1 if XIII_1_e == 2
replace obj_cal_bie = 1 if xiii_1_e == 1
replace obj_cal_bie = 1 if xiii_1_e == 2

label var obj_cal_bie "1 si la importancia de mejorar la calidad de los bie o serv es media o alta "
 
*Objetivo - Reducir los costos de producci�n

gen obj_red_cp=0

replace obj_red_cp =1 if XIII_1_h== 1
replace obj_red_cp =1 if XIII_1_h== 2
replace obj_red_cp =1 if xiii_1_h== 1
replace obj_red_cp =1 if xiii_1_h== 2

label var obj_red_cp "1 si la importancia de reducir los cost. prod es media o alta "



*Vaeriables de total de trabajadores**

gen tot_PHD=.
replace tot_PHD=XIV_1_a  if a�o==2011
replace tot_PHD=xiv_1_a  if a�o==2014

gen tot_maestria=.
replace tot_maestria=XIV_1_b   if a�o==2011
replace tot_maestria=xiv_1_b  if a�o==2014

gen tot_especialista=.
replace tot_especialista=XIV_1_c   if a�o==2011
replace tot_especialista=xiv_1_c  if a�o==2014

gen tot_tercer_nivel=.
replace tot_tercer_nivel=XIV_1_d   if a�o==2011
replace tot_tercer_nivel=xiv_1_d  if a�o==2014

gen tot_tecnico=.
replace tot_tecnico=XIV_1_e   if a�o==2011
replace tot_tecnico=xiv_1_e  if a�o==2014

gen tot_secundaria=.
replace tot_secundaria=XIV_1_f   if a�o==2011
replace tot_secundaria=xiv_1_f  if a�o==2014

gen tot_primaria=.
replace tot_primaria=XIV_1_g   if a�o==2011
replace tot_primaria=xiv_1_g  if a�o==2014
 
label var tot_PHD "total de trabajadores con calificaci�n PHD"
label var tot_maestria "total de trabajadores con calificaci�n maestr�a"
label var tot_especialista "total de trabajadores con calificaci�n especialista"
label var tot_tercer_nivel "total de trabajadores con calificaci�n tercer_nivel"
label var tot_tecnico "total de trabajadores con calificaci�n t�cnico o t. superior"
label var tot_secundaria "total de trabajadores con calificaci�n secundaria"
label var tot_primaria "total de trabajadores con calificaci�n primaria"

********% NIVEL EDUCATIVO DE LOS EMPLEADOS********

replace num_emple=1 if num_emple==0

gen PHD=.
replace PHD=XIV_1_a/num_emple  if a�o==2011
replace PHD=xiv_1_a/num_emple  if a�o==2014

gen maestria=.
replace maestria=XIV_1_b/num_emple   if a�o==2011
replace maestria=xiv_1_b/num_emple  if a�o==2014

gen especialista=.
replace especialista=XIV_1_c/num_emple   if a�o==2011
replace especialista=xiv_1_c/num_emple  if a�o==2014

gen tercer_nivel=.
replace tercer_nivel=XIV_1_d/num_emple   if a�o==2011
replace tercer_nivel=xiv_1_d/num_emple  if a�o==2014

gen tecnico=.
replace tecnico=XIV_1_e/num_emple   if a�o==2011
replace tecnico=xiv_1_e/num_emple  if a�o==2014

gen secundaria=.
replace secundaria=XIV_1_f/num_emple   if a�o==2011
replace secundaria=xiv_1_f/num_emple  if a�o==2014

gen primaria=.
replace primaria=XIV_1_g/num_emple   if a�o==2011
replace primaria=xiv_1_g/num_emple  if a�o==2014
 
label var PHD "% de trabajadores con calificaci�n PHD"
label var maestria "% de trabajadores con calificaci�n maestr�a"
label var especialista "% de trabajadores con calificaci�n especialista"
label var tercer_nivel "% de trabajadores con calificaci�n tercer_nivel"
label var tecnico "% de trabajadores con calificaci�n t�cnico o t. superior"
label var secundaria "% de trabajadores con calificaci�n secundaria"
label var primaria "% de trabajadores con calificaci�n primaria"

 
 *****LOGARITMO DE LA EDUCACION******

gen l_PHD=.
replace l_PHD=log(XIV_1_a+1)  if a�o==2011
replace l_PHD=log(xiv_1_a+1) if a�o==2014

gen l_maestria=.
replace l_maestria=log(XIV_1_b+1)   if a�o==2011
replace l_maestria=log(xiv_1_b+1)  if a�o==2014

gen l_especialista=.
replace l_especialista=log(XIV_1_c+1)  if a�o==2011
replace l_especialista=log(xiv_1_c+1)  if a�o==2014

gen l_tercer_nivel=.
replace l_tercer_nivel=log(XIV_1_d+1)   if a�o==2011
replace l_tercer_nivel=log(xiv_1_d+1)  if a�o==2014

gen l_tecnico=.
replace l_tecnico=log(XIV_1_e+1)  if a�o==2011
replace l_tecnico=log(xiv_1_e+1)  if a�o==2014

gen l_secundaria=.
replace l_secundaria=log(XIV_1_f+1)   if a�o==2011
replace l_secundaria=log(xiv_1_f+1)  if a�o==2014

gen l_primaria=.
replace l_primaria=log(XIV_1_g+1)   if a�o==2011
replace l_primaria=log(xiv_1_g+1)  if a�o==2014

label var l_PHD "logaritmo del total trabajadores con calificaci�n PHD"
label var l_maestria "logaritmo del total trabajadores con calificaci�n maestr�a"
label var l_especialista "logaritmo del total trabajadores con calificaci�n especialista"
label var l_tercer_nivel "logaritmo del total trabajadores con calificaci�n tercer_nivel"
label var l_tecnico "logaritmo del total trabajadores con calificaci�n t�cnico o t. superior"
label var l_secundaria "logaritmo del total trabajadores con calificaci�n secundaria"
label var l_primaria "logaritmo del total trabajadores con calificaci�n primaria"

 
**Postgrado (PHD y Maestr�a)

gen postgrado = 0
replace postgrado = XIV_1_a + XIV_1_b if a�o == 2011
replace postgrado = xiv_1_a + xiv_1_b if a�o == 2014

label var postgrado "total de trabajadores con calificaci�n postgrado"

**Logarirmo de postgrado

gen l_postgrado=0

replace l_postgrado=log(postgrado+1)

label var l_postgrado "logaritmo del total trabajadores con calificaci�n postgrado"

**Porcentaje postgrado

gen porc_postgrado=0
replace porc_postgrado=postgrado/num_emple 

label var l_postgrado "% de trabajadores con calificaci�n postgrado"



****Intensidades

*1. Intensidad en I+D

	*int_ga_IyD

*2. Intensidad de gastos en innovacione de productos y proceso

gen int_ga_inn_pyp=0

replace int_ga_inn_pyp = ga_inn_pyp/(vent_pn_eym + vent_pn_e + vent_p_na)

label var int_ga_inn_pyp "int. gastos en innovacion de prod y procesos (gastos en inn prod y proce / ingresos por ventas"
 

******AGRUPACION DE OBJETIVOS DE COLABORACION******

*CAPACITACION Y ASISTENCIA TECNICA

gen CyAT=0
replace CyAT=1 if Capac==1
replace CyAT=1 if AsTec==1

label var CyAT "1 si el objetivo es capacitacion o asistencia tecnica"

*I+D, INGENIERIA Y DISE�O Y PRUEBA PRODUCTOS

gen inno=0
replace inno=1 if IyD==1
replace inno=1 if IngDise�o==1
replace inno=1 if PruProd==1

label var inno "1 si el objetivo es I+D, ingenieria y dise�o o prueba de productos"

*INFORMACION

	*variable Info

*FINANCIAMIENTO

	*variable Financia


	
*******VARIABLES DE COLABORACION AMPLIADAS POR OBJETIVOS*********

*Clientes

gen clie_CyAT=0
replace clie_CyAT=1 if col_clientes==1 & CyAT==1

gen clie_inno=0
replace clie_inno=1 if col_clientes==1 & inno==1

gen clie_info=0
replace clie_info=1 if col_clientes==1 & Info==1

gen clie_fin=0
replace clie_fin=1 if col_clientes==1 & Financia==1

*Proveedores

gen prove_CyAT=0
replace prove_CyAT=1 if col_prove==1 & CyAT==1

gen prove_inno=0
replace prove_inno=1 if col_prove==1 & inno==1

gen prove_info=0
replace prove_info=1 if col_prove==1 & Info==1

gen prove_fin=0
replace prove_fin=1 if col_prove==1 & Financia==1


*Cooperacion horizontal

gen horiz_CyAT=0
replace horiz_CyAT=1 if coop_horiz==1 & CyAT==1

gen horiz_inno=0
replace horiz_inno=1 if coop_horiz==1 & inno==1

gen horiz_info=0
replace horiz_info=1 if coop_horiz==1 & Info==1

gen horiz_fin=0
replace horiz_fin=1 if coop_horiz==1 & Financia==1

*Cooperacion con entidades de conocimiento

gen conoc_CyAT=0
replace conoc_CyAT=1 if coop_conoc==1 & CyAT==1

gen conoc_inno=0
replace conoc_inno=1 if coop_conoc==1 & inno==1

gen conoc_info=0
replace conoc_info=1 if coop_conoc==1 & Info==1

gen conoc_fin=0
replace conoc_fin=1 if coop_conoc==1 & Financia==1


label var clie_CyAT "1 si colabora con clientes con obj. de Capac o Asis Tecn"
label var clie_inno "1 si colabora con clientes con obj. de I+D, ing y dise�o o prueba de prod"
label var clie_info "1 si colabora con clientes con obj. de Informacion"
label var clie_fin "1 si colabora con clientes con obj. de Financiamiento"
label var prove_CyAT "1 si colabora con proveed con obj. de Capac o Asis Tecn"
label var prove_inno "1 si colabora con proveed con obj. de I+D, ing y dise�o o prueba de prod"
label var prove_info "1 si colabora con proveed con obj. de Informacion"
label var prove_fin "1 si colabora con proveed con obj. de Financiamiento"
label var horiz_CyAT "1 si realiza coop horiz con obj. de Capac o Asis Tecn"
label var horiz_inno "1 si realiza coop horiz con obj. de I+D, ingenieria y dise�o o prueba de productos"
label var horiz_info "1 si realiza coop horiz con obj. de Informacion"
label var horiz_fin "1 si realiza coop horiz con obj. de Financiamiento"
label var conoc_CyAT "1 si colabora con ent de conoc con obj. de Capac o Asis Tecn"
label var conoc_inno "1 si colabora con ent de conoc con obj. de I+D, ingenieria y dise�o o prueba de productos"
label var conoc_info "1 si colabora con ent de conoc con obj. de Informacion"
label var conoc_fin "1 si colabora con ent de conoc con obj. de Financiamiento"
	
	

	**AJUSTE 3 	- redefinir fin de la cooperaci�n (a) innovaci�n  , b) lo demas (CyAT, informaci�n y financiamiento)
	
gen otro_obj=0
replace otro_obj=1 if CyAT==1
replace otro_obj=1 if Info==1
replace otro_obj=1 if Financia==1

label var otro_obj "1 si el objetivo es diferente a innovaci�n (I+D, ingenieria y dise�o o prueba de productos)"	

*AJUSTE 4 - REDEFINIR EL CRUCE DE COOPERANTES Y FIN DE COOPERACION

*******VARIABLES DE COLABORACION AMPLIADAS POR OBJETIVOS*********

*COOPERANTES DE PROCESO

gen proc_otro_obj=0
replace proc_otro_obj=1 if coop_proc==1 & otro_obj==1

gen proc_inno=0
replace proc_inno=1 if coop_proc==1 & inno==1

label var proc_otro_obj "1 si realiza cooperacion en proceso con objetivo distinto a innovacion"
label var proc_inno "1 si realiza cooperacion en proceso con objetivo de innovacion (I+D, ing y dise�o o prueba de prod)"


*COOPERANTES DE CONOCIMIENTO 
 
gen conoc_otro_obj=0
replace conoc_otro_obj=1 if coop_conoc==1 & otro_obj==1

*variable conoc_inno

label var conoc_otro_obj "1 si realiza cooperacion con ent. conocimiento con objetivo distinto a innovacion"
*label var conoc_inn "1 si realiza cooperacion con ent. conocimiento con objetivo de innovacion (I+D, ing y dise�o o prueba de prod)"
 

 
 
 
***VARIABLES COOPERA CON OBJETIVO DE INNOVACION Y COOPERA CON OTRO OBJETIVO (CATEGORIA BASE ES NO COOPERA)

gen coop_inno=0
replace coop_inno=1 if coopera==1 & inno==1
label var coop_inno "1 si coopera con objetivo de innovacion"


gen coop_otr_obj=0
replace coop_otr_obj=1 if coopera==1 & otro_obj==1
label var coop_otr_obj "1 si coopera con un objetivo distinto a innovacion"




***AJUSTE 30/6/2021 CATEGORIAS CON OBST�CULOS DE PERSONAL CALIFICADO EN EMPRESA Y PA�S

gen pc_categ1=0
replace pc_categ1=1 if obs_pers_cal==1 & obst_pers_pais==0

gen pc_categ2=0
replace pc_categ2=1 if obs_pers_cal==0 & obst_pers_pais==0

gen pc_categ3=0
replace pc_categ3=1 if obs_pers_cal==0 & obst_pers_pais==1 

gen pc_categ4=0 
replace pc_categ4=1 if obs_pers_cal==1 & obst_pers_pais==1


label var pc_categ1 "1 si el obstaculo de personal cal. en la empresa es alto y en el pa�s es bajo"
label var pc_categ2 "1 si el obstaculo de personal cal. en la empresa es bajo y en el pa�s es bajo"
label var pc_categ3 "1 si el obstaculo de personal cal. en la empresa es bajo y en el pa�s es alto"
label var pc_categ4 "1 si el obstaculo de personal cal. en la empresa es alto y en el pa�s es alto"


*Cruce obst�culo socio con Coopera en general

gen s_coop_categ1=0
replace s_coop_categ1=1 if obst_socio==1 & coopera==0

gen s_coop_categ2=0
replace s_coop_categ2=1 if obst_socio==0 & coopera==0

gen s_coop_categ3=0
replace s_coop_categ3=1 if obst_socio==0 & coopera==1

gen s_coop_categ4=0
replace s_coop_categ4=1 if obst_socio==1 & coopera==1

label var s_coop_categ1 "1 si el obstaculo de conseguir socio es alto y la empresa no coopera"
label var s_coop_categ2 "1 si el obstaculo de conseguir socio es bajo y la empresa no coopera"
label var s_coop_categ3 "1 si el obstaculo de conseguir socio es bajo y la empresa coopera"
label var s_coop_categ4 "1 si el obstaculo de conseguir socio es alto y la empresa coopera"



*Cruce obst�culo fondos y financiamiento publico


gen f_fp_categ1=0
replace f_fp_categ1=1 if obst_fondos==1 & fin_publ==0

gen f_fp_categ2=0
replace f_fp_categ2=1 if obst_fondos==0 & fin_publ==0

gen f_fp_categ3=0
replace f_fp_categ3=1 if obst_fondos==0 & fin_publ==1

gen f_fp_categ4=0
replace f_fp_categ4=1 if obst_fondos==1 & fin_publ==1

label var f_fp_categ1 "1 si el obstaculo de conseguir fondos es alto y la empresa no recibe fin publico"
label var f_fp_categ2 "1 si el obstaculo de conseguir fondos es bajo y la empresa no recibe fin publico"
label var f_fp_categ3 "1 si el obstaculo de conseguir fondos es bajo y la empresa recibe fin publico"
label var f_fp_categ4 "1 si el obstaculo de conseguir fondos es alto y la empresa recibe fin publico"




****ESTADISTICA DESCRIPTIVA*****

*DIFERENCIA DE MEDIAS

ttest porc_vent_pn_eym==porc_vent_pn_e  , unpaired
ttest porc_vent_pn_eym==porc_vent_p_na  , unpaired
ttest porc_vent_p_na==porc_vent_pn_e  , unpaired


tabstat porc_vent_pn_eym, statistics( mean count min max sd ) by(coop_inno)    
tabstat porc_vent_pn_eym, statistics( mean count min max sd ) by(coop_otr_obj)    

tabstat porc_vent_pn_e, statistics( mean count min max sd ) by(coop_inno)
tabstat porc_vent_pn_e, statistics( mean count min max sd ) by(coop_otr_obj)    
    
tabstat porc_vent_p_na, statistics( mean count min max sd ) by(coop_inno)    
tabstat porc_vent_p_na, statistics( mean count min max sd ) by(coop_otr_obj)    

tabstat ga_IyD, statistics( mean count min max sd ) by(coop_inno)    
tabstat ga_IyD, statistics( mean count min max sd ) by(coop_otr_obj)    



tab coop_inno coop_otr_obj



***Para diferencia de medias (art�culo de 10 hojas)***

gen inn_eym=.
gen otr_eym=.
gen no_eym=.
replace inn_eym= porc_vent_pn_eym if coop_inno==1
replace otr_eym = porc_vent_pn_eym if coop_otr_obj==1
replace no_eym = porc_vent_pn_eym if coopera == 0
gen inn_e=.
gen otr_e=.
gen no_e=.
replace inn_e = porc_vent_pn_e if coop_inno==1
replace otr_e = porc_vent_pn_e if coop_otr_obj==1
replace no_e = porc_vent_pn_e if coopera == 0
gen inn_na=.
gen otr_na=.
gen no_na=.
replace inn_na= porc_vent_p_na if coop_inno==1
replace otr_na= porc_vent_p_na if coop_otr_obj==1
replace no_na= porc_vent_p_na if coopera==0
gen inn_IyD=.
gen otr_IyD=.
gen no_IyD=.
replace inn_IyD=ga_IyD if coop_inno==1
replace otr_IyD=ga_IyD if coop_otr_obj==1
replace no_IyD=ga_IyD if coopera==0

ttest inn_eym == otr_eym  , unpaired
ttest inn_e == otr_e  , unpaired
ttest inn_na == otr_na  , unpaired
ttest inn_IyD == otr_IyD  , unpaired

ttest inn_eym == no_eym  , unpaired
ttest inn_e == no_e  , unpaired
ttest inn_na == no_na  , unpaired
ttest inn_IyD == no_IyD  , unpaired

ttest otr_eym == no_eym  , unpaired
ttest otr_e == no_e  , unpaired
ttest otr_na == no_na  , unpaired
ttest otr_IyD == no_IyD  , unpaired



***Diferencia de medias articulo - para el 5 de mayo de 2021

gen IyD_eym=.
gen IyD_e=.
gen IyD_na=.
gen IyD_gasIyD=.

replace IyD_eym = porc_vent_pn_eym if IyD==1
replace IyD_e = porc_vent_pn_e if IyD==1
replace IyD_na = porc_vent_p_na if IyD==1
replace IyD_gasIyD = ga_IyD if IyD==1

ttest IyD_eym == no_eym  , unpaired
ttest IyD_e == no_e  , unpaired
ttest IyD_na == no_na  , unpaired
ttest IyD_gasIyD == no_IyD  , unpaired

tabstat porc_vent_pn_eym, statistics( mean count min max sd ) by( IyD )
tabstat porc_vent_pn_e, statistics( mean count min max sd ) by( IyD )
tabstat porc_vent_p_na, statistics( mean count min max sd ) by( IyD )
tabstat ga_IyD, statistics( mean count min max sd ) by( IyD )

tabstat porc_vent_pn_eym, statistics( mean count min max sd ) by( conoc_CyAT )
tabstat porc_vent_pn_e, statistics( mean count min max sd ) by( conoc_CyAT )
tabstat porc_vent_p_na, statistics( mean count min max sd ) by( conoc_CyAT )
tabstat ga_IyD, statistics( mean count min max sd ) by( conoc_CyAT )

tabstat porc_vent_pn_eym, statistics( mean count min max sd ) by( IngDise�o )
tabstat porc_vent_pn_e, statistics( mean count min max sd ) by( IngDise�o )
tabstat porc_vent_p_na, statistics( mean count min max sd ) by( IngDise�o )
tabstat ga_IyD, statistics( mean count min max sd ) by( IngDise�o )

tabstat porc_vent_pn_eym, statistics( mean count min max sd ) by( PruProd   )
tabstat porc_vent_pn_e, statistics( mean count min max sd ) by( PruProd   )
tabstat porc_vent_p_na, statistics( mean count min max sd ) by( PruProd   )
tabstat ga_IyD, statistics( mean count min max sd ) by( PruProd   )

tabstat porc_vent_pn_eym, statistics( mean count min max sd ) by( conoc_info )
tabstat porc_vent_pn_e, statistics( mean count min max sd ) by( conoc_info )
tabstat porc_vent_p_na, statistics( mean count min max sd ) by( conoc_info )
tabstat ga_IyD, statistics( mean count min max sd ) by( conoc_info )

tabstat porc_vent_pn_eym, statistics( mean count min max sd ) by( conoc_fin )
tabstat porc_vent_pn_e, statistics( mean count min max sd ) by( conoc_fin )
tabstat porc_vent_p_na, statistics( mean count min max sd ) by( conoc_fin )
tabstat ga_IyD, statistics( mean count min max sd ) by( conoc_fin )




tab IyD micro
tab IyD peque�a
tab IyD mediana
tab IyD grande

tab Financia micro
tab Financia peque�a
tab Financia mediana
tab Financia grande




*Diferencai de medias Jaime

gen rad_publ=.
gen inc_publ=.
gen rad_priv=.
gen inc_priv=.
replace rad_publ = porc_vent_pn_eym if publica==1
replace inc_publ = porc_vent_pn_e if publica==1
replace rad_priv = porc_vent_pn_eym if publica==0
replace inc_priv = porc_vent_pn_e if publica==0

ttest rad_publ == rad_priv  , unpaired
ttest inc_publ == inc_priv  , unpaired
ttest IyD_eym == no_eym  , unpaired
ttest IyD_eym == no_eym  , unpaired






 
 *************************************MODELOS PARA EL ARTICULO**************************************************

**AJUSTES 30/6/2021   -     (categor�as y cruces)

tobit porc_vent_pn_e conoc_CyAT IyD IngDise�o PruProd conoc_info conoc_fin grup_empr sect_manuf micro peque�a mediana l_ga_IyD porc_postgrado tercer_nivel secundaria exp exp_2 year exporta l_ga_inn_pyp publica obst_cost_alt obst_fin_ext obst_inf_tecno pc_categ1 pc_categ2 pc_categ3 s_coop_categ4 s_coop_categ2 s_coop_categ3 f_fp_categ4 f_fp_categ2 f_fp_categ3 ,ll vce(robust)
estimates store incremental_innovdesg
margins, dydx(*) predict(ystar(0,.)) post
eststo mar_incremental_innovdesg

tobit porc_vent_pn_eym conoc_CyAT IyD IngDise�o PruProd conoc_info conoc_fin grup_empr sect_manuf micro peque�a mediana l_ga_IyD porc_postgrado tercer_nivel secundaria exp exp_2 year exporta l_ga_inn_pyp publica obst_cost_alt obst_fin_ext obst_inf_tecno pc_categ1 pc_categ2 pc_categ3 s_coop_categ4 s_coop_categ2 s_coop_categ3 f_fp_categ4 f_fp_categ2 f_fp_categ3 ,ll vce(robust)
estimates store radical_innovdesg
margins, dydx(*) predict(ystar(0,.)) post
eststo mar_radical_innovdesg


tobit porc_vent_p_na conoc_CyAT IyD IngDise�o PruProd conoc_info conoc_fin grup_empr sect_manuf micro peque�a mediana l_ga_IyD porc_postgrado tercer_nivel secundaria exp exp_2 year exporta l_ga_inn_pyp publica obst_cost_alt obst_fin_ext obst_inf_tecno pc_categ1 pc_categ2 pc_categ3 s_coop_categ4 s_coop_categ2 s_coop_categ3 f_fp_categ4 f_fp_categ2 f_fp_categ3 ,ll vce(robust)
estimates store nomejora_innovdes
margins, dydx(*) predict(ystar(0,.)) post
eststo mar_nomejora_innovdes


outreg2 [incremental_innovdesg radical_innovdesg nomejora_innovdes] using prueba18.doc, replace 
esttab incremental_innovdesg radical_innovdesg nomejora_innovdes using  "C:\Users\USUARIO\Desktop\Nueva carpeta(2)\Tobit_modelos_con_obst.csv"
esttab mar_incremental_innovdesg mar_radical_innovdesg mar_nomejora_innovdes using "C:\Users\USUARIO\Desktop\Nueva carpeta(2)\1efectos marginales_modelos_con_obst.csv" 


 
 
 
 
 
 
 
 
 
 

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
***********************TESIS***********************








*******MODELOS FINAL CORREGIDO - EXPERIENCIA Y DESTINO DE COMPRA


reg porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp latinoamerica eeuu_can europa asia africa oceania obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria exp exp_2 year
imtest, white
*se corrige heterosedasticidad
reg porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp latinoamerica eeuu_can europa asia africa oceania obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria exp exp_2 year, vce(robust)

predict res21
histogram res21, normal
ovtest
vif
*Pasa los test de correcta especificaci�n y multicolinealidad
reg porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp latinoamerica eeuu_can europa asia africa oceania obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria exp exp_2 year, vce(robust)
eststo mcopn_eym

margins, dydx(*) post 
eststo mcomarpn_eym



tobit porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp latinoamerica eeuu_can europa asia africa oceania obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria exp exp_2 year,ll vce(robust)
eststo tobitpn_eym

esttab mcopn_eym tobitpn_eym

margins, dydx(*) predict(ystar(0,.)) post
eststo tobmarpn_eym

esttab mcomarpn_eym tobmarpn_eym
 

**Histograma predicciones
 
reg porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp latinoamerica eeuu_can europa asia africa oceania obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria exp exp_2 year, vce(robust)
predict y_hat_mco
histogram y_hat_mco 

tobit porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp latinoamerica eeuu_can europa asia africa oceania obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria exp exp_2 year,ll vce(robust)
predict y_hat_tobit, e(0,.)
histogram y_hat_tobit
 
 ***EXPORTANDO RESULTADOS**

outreg2 [mcopn_eym tobitpn_eym] using prueba.doc, replace see beta

esttab mcomarpn_eym tobmarpn_eym using "C:\Users\USUARIO\Desktop\Nueva carpeta (2)\efectos marginales21.csv" 





****Modelo con interacciones

tobit porc_vent_pn_eym coopera sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp latinoamerica eeuu_can europa asia africa oceania obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp int_coop_pc int_coop_ci int_coop_ff int_coop_pm int_coop_cb int_coop_rcp l_postgrado l_tercer_nivel l_secundaria exp exp_2 year,ll vce(robust)













**********MODELO FINAL*******

*1.1 MODELO 1.1 (y = v_inn)

reg v_inn col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year
imtest, white
*se corrige heterosedasticidad
reg v_inn col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year, vce(robust)

predict res20
histogram res20, normal
ovtest
vif
*Pasa los test de correcta especificaci�n y multicolinealidad
reg v_inn col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year, vce(robust)
eststo mcov_inn

margins, dydx(*) post 
eststo mcomarv_inn


tobit v_inn col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year,ll vce(robust)
eststo tobitv_inn

esttab mcov_inn tobitv_inn

margins, dydx(*) predict(ystar(0,.)) post
eststo tobmarv_inn

esttab mcomarv_inn tobmarv_inn


***Prueba del Modelo 1.1

gen prob_v_inn=0
replace prob_v_inn = 1 if v_inn>0

tobit v_inn col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year,ll vce(robust)

probit prob_v_inn col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year

***se divide las estimaciones del tobit sobre la desv estandar y se compara con las estimaciones del probit 


*1.2 MODELO 1.2 (y = porc_vent_pn_eym)

reg porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year
imtest, white
*se corrige heterosedasticidad
reg porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year, vce(robust)

predict res21
histogram res21, normal
ovtest
vif
*Pasa los test de correcta especificaci�n y multicolinealidad
reg porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year, vce(robust)
eststo mcopn_eym

margins, dydx(*) post 
eststo mcomarpn_eym



tobit porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year,ll vce(robust)
eststo tobitpn_eym

esttab mcopn_eym tobitpn_eym

margins, dydx(*) predict(ystar(0,.)) post
eststo tobmarpn_eym

esttab mcomarpn_eym tobmarpn_eym


***Prueba del Modelo 1.1
gen prob_porc_pn_eym=0
replace prob_porc_pn_eym=1 if porc_vent_pn_eym>0

tobit porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year,ll vce(robust)
probit prob_porc_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year 

***se divide las estimaciones del tobit sobre la desv estandar y se compara con las estimaciones del probit 

display  .108772 / .5742932
display .2111705 - .18940151
display  -.1054762/ .5742932
display  -.2125614  - -.18366263
display .0629599 /  .5742932
display   .105639 - .10963024
display  .0473598 /  .5742932
display  .0416634 - .08246624
display   .1525945 /   .5742932
display  .3212132 - .26570835
display  -.1343447 /  .5742932
display  -.2676274 - -.23393051
display   .015393 /   .5742932
display .0338118 - .02680338
display   -.0536006  / .5742932
display  -.0885951 - -.09333316
display  .1277024  /  .5742932
display  .2554623 - .22236446
display .1584009 / .5742932
display  .2392657 - .27581887
display  .0559131 / .5742932
display  .1400502 - .09735985
display .0482826 /  .5742932
display  .1097206 - .08407308
display  -.0369571 / .5742932
display  -.07021 - -.06435232
display  -.13087 / .5742932
display  -.2398966 - -.22788011

*Prueba predicciones modelo tobit

*MCO
reg porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year, vce(robust)
predict y_hat_mco
histogram y_hat_mco
*Tobit
tobit porc_vent_pn_eym col_clientes col_prove coop_horiz coop_conoc sect_manuf sect_minas micro peque�a mediana l_ga_IyD int_exp obs_pers_cal obst_ci obst2 obj_part_m obj_cal_bie obj_red_cp l_postgrado l_tercer_nivel l_secundaria year,ll vce(robust)
predict y_hat_tobit, e(0,.)
histogram y_hat_tobit


***EXPORTANDO RESULTADOS**

outreg2 [mcov_inn tobitv_inn mcopn_eym tobitpn_eym] using prueba.doc, replace see beta

esttab mcomarv_inn tobmarv_inn mcomarpn_eym tobmarpn_eym using "C:\Users\USUARIO\Desktop\Nueva carpeta (2)\efectos marginales22.csv" 
esttab mcomarv_inn tobmarv_inn mcomarpn_eym tobmarpn_eym using "C:\Users\USUARIO\Desktop\Nueva carpeta (2)\efectos marginalesaa.doc"




